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ABSTRACT 

We address the issue of electromagnetic pulsar spindown by combining our experience from the 
two limiting idealized cases which have been studied in great extent in the past: that of an aligned 
rotator where ideal MHD conditions apply, and that of a misaligned rotator in vacuum. We construct 
a spindown formula that takes into account the misalignment of the magnetic and rotation axes, and 
the magnetospheric particle acceleration gaps. We show that near the death line aligned rotators spin 
down much slower than orthogonal ones. In order to test this approach, we use a simple Monte Carlo 
method to simulate the evolution of pulsars and find a good ht to the observed pulsar distribution in 
the P — P diagram without invoking magnetic field decay. Our model may also account for individual 
pulsars spinning down with braking index n < 3, by allowing the corotating part of the magnetosphere 
to end inside the light cylinder. We discuss the role of magnetic reconnection in determining the pulsar 
braking index. We show, however, that n ~ 3 remains a good approximation for the pulsar population 
as a whole. Moreover, we predict that pulsars near the death line have braking index values n > 3, and 
that the older pulsar population has preferentially smaller magnetic inclination angles. We discuss 
possible signatures of such alignment in the existing pulsar data. 

Subject headings: MHD - Pulsars 


1. INTRODUCTION 

The current canonical pulsar paradigm is that of a 
magnetized rotating neutron star (see Mestel 1999 for 
a review). However, we feel that certain fundamental 
aspects of the paradigm still remain unclear. 

One aspect of the paradigm which we hope to elucidate 
in the present work has to do with the way the neutron 
star spins down. A spinning down neutron star with mass 
M,, radius r*, and angular velocity fl loses rotational 
kinetic energy at a rate 

L = . (I) 

Here, M* ~ I.4M0, r* ~ 10 km, and (...) = d(...)/dt. 
Energy is lost through electromagnetic torques in the 
magnetosphere, although other physical processes have 
at times also been discussed (gravitational radiation, 
wind outflow, star-disk interaction, etc.). 

To a first approximation, the stellar magnetic field may 
be considered as that of a rotating magnetic dipole. Even 
under such a simplification the general description of the 
stellar magnetosphere is a formidable three dimensional 
problem, since, in general, the magnetic and rotation 
axes do not coincide. Awaiting the development of the 
general theory, one can still derive important conclusions 
based on two idealized limiting cases which have been 
studied in great extent: the case of an aligned magnetic 
dipole rotating in an atmosphere with freely available 
electric charges (i.e. with ideal MHD conditions), and 
that of a misaligned magnetic dipole rotating in vacuum. 

The neutron star is not surrounded by vacuum, and 
one needs to take into consideration the electric fields 
that develop and the electric currents that flow in 
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the rotating charged magnetosphere (Goldreich & Ju¬ 
lian 1969). The most recent numerical calculation of the 
simplest possible case, that of the magnetosphere of an 
aligned rotator in force-free approximation (Contopou- 
los 2005, hereafter C05), yielded the following rather gen¬ 
eral result for the electromagnetic energy loss 


L 


aligned 


6c 


( 2 ) 


(the quantities H f and ipopen are defined below). In that 
picture, the magnetosphere consists of a corotating re¬ 
gion of closed fieldlines which extends up to a distance 
Vc from the rotation axis, and an open fieldline region 
with enclosed magnetic flux 

B-dS = I.23:|^, (3) 

27r Jopen 2rc 


where is the polar value of the magnetic field (Con- 
topoulos, Kazanas & Fendt 1999, hereafter CKF; Gruzi- 
nov 2005; C05; Timokhin 2005; see also Appendix A). 
The above expression is valid when Tc r*. In the 
limit Tc = r*, straightforward calculation yields V'open = 
B,rll2. 

The neutron star spins down because of the establish¬ 
ment of a large scale poloidal electric current circuit flow¬ 
ing along open held lines, and returning along the edge 
of the open held line region (see CKF for a detailed de¬ 
scription). The electric current that hows between the 
magnetic axis (characterized hy ip = 0) and the edge of 
the open held line region (characterized hy ip = ipopea) 
generates the spindown torque which leads to eq. [3 The 
quantity 11^7’ in eq.|3is the angular frequency of rotation 
of the open held lines. It is set by the electric potential 
drop that develops accross open held lines, between the 
magnetic axis and the edge of the open held line region. 
This potential in the magnetosphere is in general smaller 
than the corresponding electric potential drop on the sur¬ 
face of the star. The diherence between the two potential 
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drops is just the particle acceleration gap potential which 
develops along open magnetic field lines in the vicinity 
of the polar cap. Consequently, the angular velocity of 
open field lines will in general be different (smaller) than 
n. Models of particle acceleration and pair creation of 
rotation-powered pulsars yield values of the gap poten¬ 
tial Vgap of the order of 10^^ Volts (e.g. Hibschmann & 
Arons 2001). One can directly show (see Appendix B) 
that the above together with the simplifying assumption 
that Oj’ is uniform accross open field lines, yield 


— n fldeath , (4) 

where, 

f^death = . (5) 

V^open 


This describes the so-called pulsar ‘death’, i.e. the stop¬ 
ping of pulsar emission. As the neutron star slows down 
and Q drops below fldeath, the gap potential cannot at¬ 
tain the value required for particle acceleration and con¬ 
sequent pulse generation, and the pulsar stops generating 
radio emission. 

A misaligned dipole rotating in vacuum loses energy at 
a rate 


L 


vacuum 




sin^ 9 , 


( 6 ) 


where 6 is the misalignment angle between the magnetic 
and rotation axes. We know that in real life the neutron 
star is not surrounded by vacuum, and we may argue 
that, in analogy to the aligned case, the magnetosphere 
consists of a corotating and an open line region. We may 
thus rewrite eq. El in a more general form that expresses 
the energy loss rate of an orthogonal {6 — 90°) magnetic 
rotator as 


-'orthogonal 


6 c 


( 7 ) 


We would like to emphasize at this point that eqs. El 
& El being so similar, led most researchers to ignore the 
dependence on 9 and fliT- in estimates of stellar magnetic 
fields, and in most studies of the P — P diagram. The 
aim of the present work is to show that the dependence 
on 9 and Qp is important and should not be ignored, 
especially in old pulsars approaching their death. As we 
will see, pulsar death manifests itself in a most interesting 
way through its dependence on flp. Moreover, the 9 
dependence ‘softens’ the distribution of pulsars around 
the death line in the P — P diagram. 


2. ELECTROMAGNETIC ENERGY LOSSES 


In real life, pulsars are neither aligned nor perpendic¬ 
ular rotators. It is natural, therefore, to expect that the 
two electromagnetic energy loss terms described in eqs. El 
&□ contribute together in the total energy loss. A pos¬ 
sible combination is: 


Torthogonal sin 9 -j- Taiigned COS 9 


2f^Vope„ , . 2n , 2 

-—— (sm 9 + cos 0) 

3c il 


4cr^ 


,(sin 2 0 + Tl-^^Vos^e) . 


( 8 ) 


We stress that this combination, although inspired by 
physical limiting cases, is not a rigorous derivation. An 


exact expression will require the development of a de¬ 
tailed non-axisymmetric MHD theory of the pulsar mag¬ 
netosphere. As we will show, however, this rather simple 
general expression accounts qualitatively for the main 
contributions in the electromagnetic energy loss mecha¬ 
nism. 

A key element in the present discussion is what deter¬ 
mines '(/’open, or equivalently, what determines the equa¬ 
torial extent of the closed line region. In the context of 
force-free ideal MHD, the only natural length scale is the 
light cylinder distance ric = c/D. We could have directly 
assumed that Tc is close to Gc, i.e. that the corotating 
closed line region extends up to the light cylinder dis¬ 
tance (see CKF). But this does not answer the question 
how the corotating region follows the light cylinder as the 
star spins down and the light cylinder moves out. Obvi¬ 
ously, in order for Tc to follow ric, the corotating closed 
line region must grow through north-south poloidal mag¬ 
netic field reconnection around r,.. One may consider two 
limiting cases: 

• Reconnection is very efficient, and the extent of the 
closed line region follows closely the light cylinder, 

i.e. 

Tc ~ ric ■ (9) 

• Reconnection is very inefficient, and the closed line 
region cannot grow. In that case, 

Tc « const. ( 10 ) 

One might argue that the supernova explosion that 
led to the formation of the spinning neutron star 
blew the magnetic field open, and the magnetic 
field remained open since then. 


We are going to consider a simple parametric model 
where rc is neither equal nor proportional to ric, but 
is equal to 


rc = ric 




( 11 ) 


where r/co and Oo are the values of ric and fl at pul¬ 
sar birth repsectively, and 0 < a < 1. Equation El bas 
a different from dipolar dependence on O, i.e. a depen¬ 
dence different from oc This accounts for braking 
index values n 7 ^ 3 (we introduce here the first and sec¬ 
ond order braking indices n = and m =ft 

respectively). 

We have presented here a natural physical model that 
accounts for braking indices n 7 ^ 3, namely a corotat¬ 
ing closed line region which does not extend up to the 
light cylinder but follows the light cylinder at an increas¬ 
ing distance. The dependence of energy loss on the ex¬ 
tent of the magnetosphere and its effect on the brak¬ 
ing index was noted in the past (Blandford & Romani 
1988; Harding et al. 1999; Alvarez & Carraminana 2004; 
Gruzinov 2005b; Timokhin 2005), with Gruzinov (2005b) 
and Timokhin (2005) introducing the same evolutionary 
parametrization as in CU- The physical motivation of 
the extent of the closed zone varies between the authors 
(e.g., inertial stresses by the wind in Blandford & Ro¬ 
mani 1988 and in Harding et al. 1999, obliqueness effects 
in Gruzinov 2005b, or inability of pulsar gaps to sup¬ 
ply necessary poloidal currents in Timokhin 2005). Our 
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Fig. 1. — P — P evolutionary diagram that shows the effect 
of the misalignment angle 6. We show spindown trajectories for 
pulsars with initial period Po = 10 msec, magnetic field = 
10^^ G, and Vgap = 10^^ Volt. Initial braking index n = 3 is 
assumed (q = 0). Trajectories from bottom to top correspond to 
pulsars with increasing magnetic inclination angle from 0° to 90° 
with 15° increments. Rectangular dots indicate when pulsars reach 
P^eath turn off. The das hed line represents the theoretical 
death line as given in eq. IT^ which does not take into account 
the misalignment angle dependence introduced in eq. mi Oblique 
pulsars evolve faster through the diagram. 


view, however, that Tc < ric might be due to insuffi¬ 
ciently fast reconnection near the Y-point in the outer 
magnetosphere is original and worthy of further study. 

More importantly, though, eq. |H1 brings a new element 
to the discussion: we directly introduce pulsar death and 
magnetic-rotation axes misalignment in the electromag¬ 
netic pulsar spindown expression. Pulsar death expresses 
the inability of the magnetosphere to generate the par¬ 
ticles required in the poloidal electric currents which 
generate spindown torques in the aligned rotator case. 
They do not affect the orthogonal spindown component 
though. Beyond the death line, the misaligned neutron 
star will continue to spin down without emitting pulsar 
radiation. 


3. PERIOD EVOLUTION AND THE P- P DIAGRAM 

In order to study how pulsars spin down, we equate 
the values of L in eqs.Ql&IHlaiid thus obtain 


P = 3.3 X 10"^® 
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where from © 


Tdeath = 8.12— sec 


ioi2Gy Vio^^voit 
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Po is the initial period at pulsar birth. Note that we 
have introduced an angular dependence in the dehnition 
of the death line P(Pdeath)- As we will see below, this has 
interesting observational consequences. Note also that, 
in the limit 9 ~ 90° and a = 0, we obtain the following 
simple expression 


Aieath = 5 X 10-1® 


-f^death 
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f ^gap \ 
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which can also be written as 



(14) 


log Pdeath = 3 log Pdeath “ 17.3 -L 2 (log Vgap - 12 ) (15) 

(compare with Zhang, Harding & Muslimov 2000). 
When P Pdeath, 

( 16 ) 

n = 3-2a (17) 

m = n{n — 1 ) , (18) 

and pulsars evolve in the logarithmic P — P diagram 
along lines of constant magnetic held P oc p 2 a-i p 
approaches Pdeath, however, the evolution curves down 
away from the above straight lines and ends at a point 
p = Pheath and P = P(Pdeath) (sce hgureip. This evo¬ 
lution is similar to what one would expect if we had as¬ 
sumed magnetic held decay in the standard dipole spin- 
down model (e.g. Gonthier et al. 2002 and 2004, here¬ 
after G02 and G04 respectively). In hgure ^ we show 
only a = 0, or evolution with initial braking index n = 3. 
For n < 3 the straight portion of the lines would rotate 
counterclockwise, and would become horizontal for n = 2 
{a = 0.5). Such pulsars would also evolve faster. 

As a test of the validity of our results, we tried to re¬ 
produce the observed distribution of pulsars in the P — P 
diagram with the minimum number of assumptions pos¬ 
sible, through a simple Monte Garlo numerical experi¬ 
ment that follows G02 and G04. We integrated eq. El 
for random times t within a 10 ® year time interval (this 
mimics random pulsar birth in the galaxy). The result of 
the integration is recorded only if the star remains active 
as a pulsar, i.e. only if P(t) < Pdeath- To make our re¬ 
sults more realistic, we assumed a multipeaked lognormal 
distribution of polar magnetic helds around P* ~ lO^^G 
(eq. linG02): pe = We used 

two terms with Ai = 20.0, A 2 = 2.0, log Pi = 12.65, 
logP 2 = 13.0, (Ti = 0.6, and (T 2 = 0.4. 

We assumed a uniform random distribution of initial 
pulse periods from 10 msec to 0.2 sec. We did account for 
the observational selection effect due to finite instrument 
detectability by implementing the radio pulsar luminos¬ 
ity model of Narayan & Ostriker 1990 (eqs. 23, 24 in 
G02), and a lower detectable pulsar luminosity limit of 
12.5 mjy-kpc®. Finally, we assumed that there is no mag¬ 
netic held or misalignment angle evolution with time. In 
an effort to simplify the approach we did not consider 
the multitude of selection effects and individual survey 
detectability limits. We also did not calculate the dis¬ 
tance to individual pulsars after tracking them in the 
galactic potential as in G02 and G04, relying on the lumi¬ 
nosity distribution function. We believe, however, that 
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the gross effects of the new spindown law C2J can be 
understood even with our procedure. 

Our results are shown in figs. 2,3 & 5. In figure |21 
we show the result of our simulation for the standard 
dipolar spindown (a = 0, or equivalently n = 3, and no 
dependence on misalignment angle 0, i.e., 9 = 90° for all 
pulsars). One can directly see a strong concentration of 
pulsars near the death line, although the concentration is 
decreased by the luminosity selection effect, which favors 
younger brighter pulsars and fills up the upper left corner 
of the diagram. This concentration is understood since 
pulsars evolve slower and slower as they grow older. This 
is, however, not what is observed (figure 4). The need to 
make the death line boundary softer led some researchers 
to introduce a ’’death valley” (G02) or field decay (G02, 
G04) in their simulations. In figure|31we show that intro¬ 
ducing the 9 dependence that we have obtained in eq. |H1 
with random misalignment angles, solves this problem 
naturally without magnetic field decay. As we saw in fig¬ 
ure Q] the angle-dependent death line helps in smoothing 
the pileup. In figure^ we show the present day observed 
radio pulsar distribution from the ATNF pulsar catalog"* 
for comparison with our results. The similarity between 
figuresEl&^is quite good considering the simplifications 
of our model. We aim to reproduce the gross features of 
the observed P — P diagram for non-millisecond pulsars 
- its triangular shape and centroid location, while keep¬ 
ing the number of initial magnetic field components to 
a minimum. In the standard spindown picture it is very 
difficult to have pulsars in the lower corner of the P — P 
triangle without polluting the low P - low P region of the 
diagram with low-field pulsars, which are not observed. 
In our spindown model the lower corner is filled with 
pulsars dropping from the higher P (higher B*) region, 
rather than the ones moving on straight tracks as in the 
standard case. This contributes to the triangular shape 
of the diagram. 

We obtain a better overall fit if we constrain the mis¬ 
alignment angles to be smaller than 45 ° (figure E)). This 
reduces the horizontal width of the triangle, which means 
that more higher-held pulsars are moving on vertical 
tracks in the diagram. Also, the ridge of overdensity near 
the theoretical death line seen in hgureElis smoothed out 
in hgureEl 

Similar shapes of the diagram are found for a > 0. 
This explains individual observed pulsars with n < 3. 
However, our ht deteriorates for values of a larger than 
0.25. We thus conclude that, although individual pul¬ 
sars may well have a > 0, as a whole, a ~ 0 is a good 
approximation for the observed pulsar population. 

Furthermore, when a = 0, from eqs. mm we derive 

2 

H/(OdeathCOs2 0) - I ■ 


Around the pulsar death line SI Sldeath, this equation 
gives 


n 


2 cos^ 9 

3 + - - 

1 — cos^ 6 


( 20 ) 


which, for aligned pulsars with cos 9^1 may be much 
greater than 3. Such high braking index values may be 
observable (Johnston & Galloway 1999). 


^ Available at http://www.atnf.csiro.au/research/pulsar/psrcat/ 
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Fig. 2.— Monte Carlo experiment using the standard dipole 
model {a = 0, or equivalently n = 3, and no angle dependence). 
Pulsars are injected with a lognormal distribution of polar magnetic 
field values around = 10^^ G (as in G02). V^ap = 10^^ Volt. 
Po, the initial period at pulsar birth, is uniformly randomly dis¬ 
tributed between 10 msec and 0.2 sec. We have also introduced an 
observational selection effect due to finite instrument detectability 
limit (see text). T he d ashed line represents the theoretical death 
line as given in eq. im The strong concentration of pulsars near 
the death line is not observed in the real distribution of pulsars. 


In figure|n|we plot the distribution of cos 9 in our Monte 
Carlo simulation that was presented in figure |31 In this 
simulation pulsars are injected (born) with random mag¬ 
netic axis orientation with respect to the rotation axis as 
shown by the dashed histogram in figure El The pul¬ 
sar population at the end of the simulation, however, 
shows a different distribution of misalignment angles as 
shown by the thick solid line. The selection bias towards 
smaller misalignment angles may be attributed to the 
fact that orthogonal rotators spin down faster (on the 
average), and therefore are more likely to evolve beyond 
the death line within the observation time interval. Fi¬ 
nally, the thin solid line represents the distribution of 
misalignment angles in the pulsars observed lying below 
the theoretical death line given by eq. El (which does 
not take into account the angle dependence introduced 
in eq. cai. One sees clearly that pulsars observed near 
or below the death line should have preferentially smaller 
misalignment angles 9 than the rest of the pulsar popu¬ 
lation. This is, we believe, an important observational 
prediction of the present work. 

4. DISCUSSION 

Awaiting the development of a detailed three dimen¬ 
sional MHD theory for the rotating neutron star magne¬ 
tosphere, we can describe a few general characteristics of 
its expected structure. 

We argue that ftp characterizes the reduced magneto- 
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Fig. 3.— Same as fig. [^1 only with the misalignment random FiG. 5.— Same as figure 2, but with misalignement angles con- 

angle dependence obtained in eq. The introduction of the angle strained to be less than 45°. 
dependence results in the softening of the death line. 
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Fig. 4.— The present day observed radio pulsar distribution. 
Data from ATNF pulsar catalog. 


spheric electric potential drop between the magnetic axis 
and the edge of the open field line region. Poloidal elec¬ 
tric currents will be generated as long as electric charges 
can be produced in the magnetosphere. These charges 
are produced in the magnetospheric polar gaps. In the 


Fig. 6.— The distribution of |cos(0)| in our Monte Carlo ex¬ 
periment in figure 3. The dashed line represents the original dis¬ 
tribution injected with uniform random misalignment angle. The 
thick solid line represents the distribution of misalignment angles 
in the simulated pulsar distribution as a whole. The thin solid line 
represents the distribution of misalignment angles for the pul sars 
observed lying below the theoretical death line given by eq. im 
An interesting observational prediction of the present work is that 
pulsars lying near or below the theoretical death line should have 
preferentially smaller misalignment angles. 
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Fig. 7. — Distribution of alignment measure in the observed 
pulsar population. Each pulsar is plotted as a circle with radius 
proportional to the quantity using the pulse width W at 

50% of the peak. Data from ATNF pulsar catalog. 


axisymmetric case, can also be thought as the an¬ 
gular velocity of rotation of open magnetic field lines, 
which is in general smaller than Q. The electric current 
that flows between the magnetic axis and the edge of the 
open field line region generates the spindown torque. As 
shown in C05, this electric current and spindown torque 
are both proportional to flp- 

We argue that this picture describes also the situation 
when 9^0. As the neutron star rotates, the magnetic 
axis moves around the axis of rotation. In the rotating 
frame of the star, however, open field lines rotate at their 
own rate — flp opposite to the direction of stellar ro¬ 
tation around the magnetic axis. Observational evidence 
for this effect can be found in the well known sub-pulse 
drift phenomenon (e.g. Beskin 1997; Rankin & Wright 
2003). In the limit Oir = O, it is as if open magnetic field 
lines are anchored on the stellar surface. In the opposite 
limit rip = 0, open field lines rotate around both the 
magnetic and rotation axes returning to the same posi¬ 
tion every rotation of the star when viewed by an inertial 
observer. For finite a phase shift about the magnetic 
axis will be accumulated after every turn of the star. A 
simple analogy to the above picture might be that of 
a lawn watering system consisting of a hose rotating at 
angular velocity Id, and a sprinkler at its end rotating 
with angular velocity — SI in the rotating frame of 
the watering hose. 

Note that the orthogonal component does not require 
the establishment of a poloidal electric current circuit in 
order for it to slow down the stellar rotation. The or¬ 
thogonal component emits spiral electromagnetic waves 
(in general Alfven waves) which travel out to infinity 
through the rotating or non-rotating open field lines. 


The main conclusions of the present work are: 

1. The electromagnetic energy loss of a pulsar de¬ 
pends not only on the surface magnetic field R* 
and on the rotational angular velocity S7, but also 
on the misalignment angle 6 and on the angular 
velocity of the open magnetic field lines SIf (which 
is in general smaller than Q). 

2. The approach to pulsar death modifies the rate of 
energy loss. The death occurs when the pulsar 
slows down sufficiently so that the angular velocity 
of the open fieldlines = 0. 

3. The energy loss close to the pulsar death is smaller 
than what is given by the standard dipolar spin- 
down formula. This effect gives a good fit between 
the theoretical and observed distributions of pul¬ 
sars near the death line, without invoking a mag¬ 
netic field decay. 

4. Our model may also account for individual pulsars 
spinning down with braking index n < 3. However, 
n ~ 3 remains a good approximation for the pulsar 
population as a whole. 

5. Pulsars near the death line have braking index val¬ 
ues n > 3. Such high braking index values may be 
observable. 

6 . Pulsars near the death line may have preferentially 
smaller inclination angles. 

A preliminary look at the ATNF pulsar catalog data 
suggests that the last point may have some observational 
support. One possible measure of the inclination angle 
of a pulsar is its fractional pulse width, or the ratio of 
the width of the pulse to the period of the pulsar. If 
the radio beam size is independent of inclination, then 
we would expect pulsars with smaller inclinations to be 
seen for a larger fraction of the period than pulsars with 
large inclinations. In order to test this hypothesis we 
took the available data for the pulse width at 50% of the 
pulse peak from the ATNF catalog (1375 pulsars with 
W 50 % 7 ^ 0). In order to be able to compare the frac¬ 
tional pulse width for pulsars of different periods, we 
have to correct for the intrinsic size of the pulsar beam. 
If the beam roughly follows the angular size of the open 
fieldlines, then the beam size falls as with increas¬ 

ing period (we are assuming that the last closed field 
line extends out to the light cylinder). Therefore, the 
quantity that should relate to the degree of alignment 
is Fliiign = {W X P^/‘^)IP = WP~^^‘^, where W is the 
measured pulse width. In figure Owe plot the observed 
pulsars as circles with the radius of the circle linearly pro¬ 
portional to Paiign- A visual inspection of the plot shows 
that there is an excess of larger pulse fractions for older 
pulsars, and in particular for pulsars near the right edge 
of the P—P triangle. This region is near the pulsar death 
line (dashed line), and therefore, this feature is particu¬ 
larly interesting. Obviously, a much better analysis needs 
to be performed, but these results are quite encouraging. 
An overabandunce of pulsars with large pulse fractions 
near the death line in a smaller pulsar sample had also 
been interpreted by Lyne and Manchester (1988) as an 
indication of alignment. 
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In an effort to simplify the fits we have assumed in 
this paper that the misalignment angle stays the same 
during the evolution, and there is no field decay. The re¬ 
ality may of course be more complicated, and there could 
be some amount of field decay, and potential alignment 
(and consequent precession) during the lifetime of a pul¬ 
sar, due to electromagnetic torques on the star. These 
effects would introduce extra degrees of freedom to fit¬ 
ting the pulsar distribution. We hope, however, that 
the physically motivated spindown law of the type in¬ 
troduced in this work would find its way into detailed 
population synthesis models. 
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APPENDIX 

APPENDIX A: CALCULATION OF OPEN MAGNETIC FLUX 

We present here expressions for '0open(^c) obtained numerically for various axisymmetric magnetospheric configura¬ 
tions using the numerical code developed in C05: 


• Standard magnetostatic dipole: 


il) 


open — 


B^rl 

2rc 


(Al) 


Note that in this limit, V'open(?'c) is defined as the amount of flux that crosses the equator beyond equatorial 
distance Tc- All field lines cross the equator perpendicularly. 


• Magnetostatic dipole with open field lines beyond re- 

V^open = 1-20^ . (A2) 

Note that this implies the existence of an equatorial azimuthal current sheet beyond rc- The magnetosphere is 
otherwise stress-free. The last open field lines crosses the equator at an angle of about 70° towards the axis (‘Y’ 
point). 

• Rotating (relativistic) dipole with open field lines beyond Tc, and smooth crossing of the light cylinder: 

V'open = 1-23:^ . (A3) 

ZT c 

This result is valid for Tc < ric- The solution requires the presence of a) a poloidal electric current along 

open field lines, b) a poloidal current sheet along the separatrix between open and closed field lines, and c) an 
equatorial azimuthal current sheet beyond rc- The last open field line crosses the equator at an angle of about 
70° towards the axis. Note the interesting similarity with the latter non-relativistic case (eq. roi . 


APPENDIX B: CALCULATION OF 


We derive here an approximate expression for flp., assuming flp = const. We perform the calculation in the aligned 
(axisymmetric) case 0 = 0. It has been shown in COS that ^Ip is related to the magnetospheric potential drop Vp 
between the magnetic axis (characterized by ')/' = 0) and the edge of the open field line region (characterized by 
V’ = V'open), namely 

/*'?/'open /•'j/'open oQe 7 1 /*’0open 

Vf= E • ds = / —Bpds = - 

J 'Ip — Q J -0—0 ^ ^ J 0 — 0 




(BI) 


c 







where s is the distance from the rotation axis; E is the electric field; Bp is the poloidal (meridional) component of the 
magnetic field. Vf is in general different from the corresponding stellar potential drop, namely 


K = i 

c , 


nd'0 = 


^^open 




The difference 


Kan =V,-Vf = 


(n — 


open 


(B2) 


(B3) 


is just the particle acceleration gap potential which develops along the magnetic field near the footpoint of open held 
lines (e.g. Beskin 1997). Therefore, 

rip = f^death ; (^^4) 


where 


f^death — 




gap 


'Ip, 


(B5) 


open 


This expression is valid as long as > fldeath- When Q < fldeath, = 0. 



